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ABSTRACT 


The purpose of this investigation was to use narrow band Ross filter 
measurements of the Fe XXV line flux around 0.185 nm wavelength and 
simultaneous broad band measurements during a solar flare to determine 
the relationship between the solar coronal abundances of iron and helium, 

A second purpose was to determine the Fe XXV ion population as a function 
of time during a flare. Broad band measurements from the UCSD propor- 
tional counter and Ross filter measurements from the Goddard Space Flight 
Center instrument, both on OSO-7, were taken of a flare at UT 0300 on 
27 October, 1972. The UCSD measurements gave the temperature and emis- 
sion measure as a function of time and allowed line and continuum to be 
separated in the Goddard measurements. The data were analyzed under two 
separate assumptions: that the electron density was high enough that a single 
temperature could characterize the continuum spectrum and the ionization 
equilibrium, and that the electron density was low so that the ion populations 
trailed the electron temperature in time. Using the imaging property of the 
Goddard experiment, we found that the density was at least 5 X 10 , and that 
the high density assumption was valid. Fitting the observed line strength to 
the measured temperature and exnission measure we found the iron abundance 

_5 

to be 1. 1 X 10 for a helium abundance of 0, 2, relative to hydrogen. This 

result can be made consistent with previously determined iron abundances 

"• 5 

around 2 X 10 but not witn much higher values. An absolute determination 
of the helium abundance by this method awaits the analysis of many good iron 
spectra. 
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1. Introduction 


The purpose of this investigation is to measure the iron line to continuum 
ratio in the wavelength range around 0.185 nm, In the past the line emission has 
been measured using crystal spectrometers (Neupert, et al . , 1967), while 
continuum measurements have been made using broad band instruments such 
as proportional counters (McKenzie, et al . , 1973), Because crystal spectro- 
meters have low efficiency, and proportional counters experience saturation 
and pulse pileup difficulties at high counting rates, the line to the continuum 
ratio has not been measured previously. This investigation uses data from the 
Goddard Space Flight Center iron line Ross balanced filter pair and the 
University of California at San Diego proportional counter, both on OSO-7. 

Once the line to continuum ratio has been determined, we use theoretical 
calculations of the line strength to estimate the relative coronal abundances 
of iron and helium and to examine the time dependent population of the helium- 
like ion, Fe XXV. 

Table I presents the thresholds for the channels of the UCSD proportional 
counter. These values are discussed in a later section. The iron line emission 
falls primarily in channel 4. We therefore fit the data in channels 3 and 5-8 to 
a model spectrum to obtain values of the flare electron temperature, T , and 
emission measure (EM), defined as: 

EM = /" nj dV, (1.1) 

Jy 

where n g is the flare electron density, and the integral is over the flare 
volume. These determinations allow the continuum at 0. 185 nm to be estimii,* 

The GSFC Ross balanced filter pair measures the flux in a narrow, well- 
defined energy range. Using the UCSD continuum determination we separate 
the line flux from the continuum within the Ross filter band pass. This gives 
us the line to continuum ratio. Using theoretical calculations of the line and 
continuum fluxeB we are then able to estimate the relationship between the 
solar iron and helium abundances. This is possible because the contribution 
of helium to the continuum flux is second in importance only to that of hydrogen. 
Again using theoretical models, the population of Fe XXV as a function of time 
is estimated from the line flux. 
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Table I 


UCSD Energy Channel Thresholds 


Channel Number 


Prelaunch Threshold at 

Threshold (nm) 10^ c/sec (nm) 


Threshold at 
2 x IQ 4 c/sec (nm) 



2 . 


Instrumentation 


I f 


Figure 1, a plot of the mass absorption coefficients of iron and manganese 
as a function of wavelength may be used in discussing the operational principle 
of the Ross balanced filters. The mass absorption coefficient for each material 
jumps dramatically at its K shell absorption edge wavelength. The transmission 
of, for example, a thin iron filter is given by 




Fe 


= e 


Fe 




< 2 . 1 ) 


where Pp e is the mass absorption coefficient in ern /g and Xis the filter thickness 
ingm/cm^, According to Figure 1, there exist values of and Xj^ n such 
that the transmissions of iron and manganese filters are equal except in the 
region between the two elements' absorption edges (0. 1?4 - 0. 190 nm). A Ross 
filter system consists of a pair of such balanced filters in front of identical 
detectors. The difference in the counting rates of the two counters is then a 
measure of the X-ray flux in the narrow passband, 0, 174 - 0. 190 nm. In 
addition to the narrow band, the Ross filter system has a second advantage. 

The effective detector passband is precisely set by the filter absorption edge 
energies, and is therefore not subject to electronics drift uncertainties. The 
electronic threshold needs only to be set at a low enough energy that all detected 
photons of wavelength shorter than 0. 19 nm will be counted. 


The Goddard Ross filter system consists of the balanced filters with their 
identical proportional counters located in back of a multigrid collimator. The 
properties of the system were measured before launch or are easily calculable. 
The collimator is continuously rastered over a 5 arc minute x 5 arc minute 
region in the mode employed for this measurement. Data are read out 368 times 
during the raster which takes 61,44 seconds to complete. 


The system response may be computed as the product of the efficiencies 
of the various components: 


2 

R = Ae e _ T t ,, cm -sec , (2.2) 

c F c eff 

2 

where A is the window area of one counter (2.2 cm , taking into account a 
preflight adjustment to bring the system into balance), e c is the calculated 
counter efficiency, is the filter transmission difference in the passband 
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(0.62), T c is the collimator transmission, and t e ^ isi the effective average 

exposure time for each point in the raster. The counter efficiency is the 

absorption of a 1. 05 cm path length of 90% xenon multiplied by the transmission 

2 

of the 0. 012 gm/cm beryllium window. This comes out to be 0. 912. The 
on-axis collimator transmission is 0, 12, and the collimator field of view is 
15.5 arc seconds, FWHM (Sterk, et al , , 1972). Since the raster line spacing 
is 20 arc seconds and the collimator FWHM is 15. 5 arc seconds, the collimator 
response must be multiplied by a factor to account for the fact that all points are 
not swept over by the same effective area. Assuming a triangular response, 
which characterizes a multigrid collimator of the type used, this factor is 0, 775. 
During a raster, each point is observed while the collimator sweeps over 31 
arc seconds, with an average of one-half the on-axis transmission during this 
period (again assuming a triangular response'. Thus 


"cff 


^ (31 arc sec /scan rate) 


(2.3) 


The scan rate is 78. 125 arc sec/sec. Then, equation (2,2) gives R - .0229 

2 - 1 
cm -sec. The conversion factor from raster counts to flux is R = 43.6 photons/ 

, 2 

(cm -sec) per raster count. 

The UCSD detector with which we are concerned is a proportional counter of 

2. 54 cm diameter filled with a one-atmosphere 90% Xe, 10% CO, gas mixture, 

2 * 2 

The counter has a 0. 115 gm/cm beryllium window of area 0. 32 cm , The 
instrument gain and resolution were measured prior to launch. Its behavior 
at high counting rates was also observed in the laboratory. The wavelength 
threshold of each channel is given in Table I, In addition to the prelaunch nominal 

4 4-i 

values we have noted the effective values at counting rates of 10 and 2x10 sec 
at the time of the observations discussed here, in the fall of 1972. These values 
will be discussed below. With the exception of channel 4, channels 3-8 are 
dominated by continuum, as opposed to line, flux. 

The response of the UCSD detector to model thermal X-ray soectra was 
determined for temperatures in the range 1-3 x 10 K. Model spectra which 
included the continuum processes of bremsstrahlung, radiative recombination, 
and two-photon decay from metastable levels in hydrogen- and helium-like ions 
were folded through the instrument efficiency and resolution functions. It is 
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convenient to use these data in the form of ratios of channel N counts to channel 3 
counts, with N>3. Channels 1 and 2 have important line contributions and 
efficiencies that change rapidly as a function of energy, so the data in these 
channels were not used in this analysis program. Sets of curves of the ratios 
were generated for counter gains equal to the nominal value divided by a 
parameter, k. Figure 2 shows the set for k = 1.06 

Foi the analysis undertaken here it was necessary to include the instrument 

response to iron line emission around 0. 185 nm. Furthermore the data included 

4 >1 

samples with detector integral counting rates in the range 1-2 x 10 sec , In 

55 

prelaunch measurements, using a Fe radioactive source (0,21 nm), we observed 

4 -1 

a gain decrease of about 5% at a total counting rate of 2 x 10 sec . Culhane, 

et al, (1967) observed a similar decrease in their OSO-4 instrument. The 

effect is attributed to a build up of positive ions in the counter gas. The resolution 

is not impaired, and we use the prelaunch values for our analysis. Also, since 

4 -1 

the instrument pulse pair resolving time is 3.2 usee, at a rate of 2 x 10 sec 
approximately 6.2% of the pulses will "pile-up, " adding their amplitudes to those 
of preceding pulses. 

To take into account the effectB described in the previous paragraph, a 

Monte Carlo program was written. The low counting rate model spectra, 

broadened and folded through the instrument efficiency were used as input. 

The pulse shape and the peak detecting operation of the electronics were properly 

taken into account. The iron line was broadened according to the instrument 

resolution and was included in the input spectrum. The output spectrum was 

- 17 

displayed in energy bins 3. 2 x 10 joules (0.2 keV) wide. These data were 
gathered into energy channels assuming a gain shift by a small factor, k, as 
described above, with an additional 5% decrease at 2 x 10^ 2ec - ^. Comparing with 
actual flare spectra resulted in good fits for k = 1. 05, A five percent gain 
deterioration after more than a year's continuous operation is not unreasonable. 

Figure 3 shows the Monte Carlo program results for an integral counting 

4-1 3-1 

rate of 2 x 10 sec , including a line counting rate of 1 x 10 sec , Other 

runs varied both counting rates. A few remarks comparing Figure 2 and Figure 3 

would be appropriate. In Figure 3, the channel 4 curve is higher than that for 

Figure 2, reflecting the inclusion of line emission. For the other channels. 

Figure 3 shows enhanced counting rates, particularly at low temperatures 


- 6 - 


and high energies, where the thermal spectrum is steepest. This indicates 
that pileup is the major effect, since a decrease in gain would lower the counting 
rates in the higher channels. The relatively small number of pileup events form 
a significant fraction of the counts where the "real" counting rate is small. The 
error bars in Figure 3 are £ 1 a. 
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3, Calculation ot Iron Line Emission 


2 1 1 

The iron line flux around 0. 185 nm consists of the Is S - ls2p P 

2 1 3 2 1 3 

resonance line, the forbidden Is S • ls2p S lino, the Is S - ls2p P 

intercombination lines, all from Fe XXV, and satellite lines from Fo XXIV 

cmlvted following inner shell excitation or dielectronic recombination. 


The resonance line emission from a coronal plasma is discussed by 

2 1 „ , 1 , 


Walker (1972). For the Is S - ls2p P transition (hereinafter called 
the "resonance" line) of Fe XXV, the energy flux at the Earth is 

f 


$ 


R~ 


tEM)A Fe a 24 5.°l X l°- 44 ( Z 2 Ti) 2 . B r 

' Lyo 


-E R /kT 


/~T 


/ -2 - 1 . 
(erg-cm -sec ). 


(3. 1) 


In equation (3. 1): 


EM - 



z = 


{z 




emission measure, 
iron abundance, 

fractional population of 24 times ionized Fe XXV 
(from Jordan, 1970), 

the ionic charge, 

the temperature averaged collision strength for the 
Lyman alpha transition in Fe XXVI, 

the resonance line oscillator strength, 

the Lyman alpha oscillator strength, 

the branching ratio of the resonance line, 

the excitation e nergy of the resonance line, 

the electron temperature (K), and 

Boltzmann's constant. 


Note that equation (3. 1) calculates the resonance line collision strength from 
that for the Lyman alpha line. If the electron density is low the ion populations 
may be characterized by a temperature, T , which differs from T. In such 
cases, equation (*. 1) may be used with a^ equal to a 2 ^(T z ). 

Following dielectronic recombination the lithium- like ion, Fe XXIV, may 

2 

be in a state ls2p n l. If it undergoes the transition, Is ni • ls2pni rather than 
aut< onizing, the emitted photon has a wavelength nearly the same as that of 
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the resonance line. For n ** 3, the "satellite lines" are generally unresolvablc 
from tlu* resonance line. The unresolved dielcctronic recombination flux may 
he related to the resonance line flux by 


UI> 


Qf$, 


R 


(3.2) 


where 


0 ' - 


5. 375 x 10 v 


(3. 3) 


2vt 


for Fe XXV. Since (z“U ). 1.6, the unresolved dielcctronic recombination 

Lyo 

satellite lines make up 16-35% of the measured resonance line flux. Equation 
(3. 3), based upon an extrapolation of the work of Shore (1969) to high Z, probably 
overestimates a. Walker (personal communication) has rccaculated a to arrive 
at (for Fe XXV) 


a 


2. 19 x 10° 


( 3. 4) 


We performed the ana’vsis using equation (3.3), however. The resonance line 

flux one would moa r. $ , is the sum of and S’ 

s K UD 

The forbidden and intercombination line emission has been discussed by 
Gabriel and Jordan (1969, 1973) and Blumenthal, et al. (1972). The sum of 
these lines is related to the resonance line by an experimentally determined 
parameter, G: 


$ I+F “ G V 


(3.4) 


Equation (3.4) uses $ since G is an experimental parameter. G has been 

s 

determined for low Z elements. Rugge and Walker (1971) found G " 1.1. 
We use this value in the absence of a measurement for iron, 

Gabriel (1972) and Bhalla et al. , (1975) have calculated the resolved 
satellite line fluxes (n = 2). For an individual line excited by dielcctronic 
recombination, 


- 9 - 


i if t 



(E r - E D )/kT 

O X 


». A r A . 

I*?**? 


Here E D and E.^ arc the resonance and satellite line energies, g is the 
R U 8 

satellite level statistical weight, and A and A are the radiative decay 

r a 

and autoionization transition probabilities. 

Gabriel rewrites (3. 5) as 



separating out the temperature dependence. The second factor is the ratio 
evaluated at T , the temperature at which resonance line emission is a 
maximum. For our broad band measurements we sum (3.6) over the satellite 
lines, all of which lie in the Ross filte.* pass band. According to B'halla et al. 
(1975), 


* d <T) = * b (T) [o. 2963 F d (T)] . 

Fjj(T) is tabulated by Bhalla et al. and shown V Figure 4. 

Because the inner shell excitation line strengths depend on the relative 
populations of Fe XXIV and Fe XXV, they are related to the resonance line 
strength through T . According to Gabriel, 

55 

*1 - a 23 _1_ „ e !l 

V a 24 1+ » f R ' 

where 6 is the probability that the Fe XXIV ion is in the proper initial state 
to be excited and emit the satellite line in question and the f 1 s are effective 
oscillator strengths. As in (3.7) we write 

$j{T,T z ) * * s (T,T z ) jo.0808 F I (T z )j . 

Fj(T z ) is plotted in Figure 4, 
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i l ' r ‘\ 


(3.5) 


(3.6) 


(3, 7) 


(3.8) 


(3. 9) 


4. Data Analysis 

Three events occurring on 27 October, 29 October, and 3 November, 

1972, were selected by Dr. Roger J. Thomas of Goddard for possible 
analysis. All three events occurred in McMath Plage Region 12094. The 
3 November burst had UCSD counting rates of 5 x 10^ -- too high for reliable 
spectral analysis. The 27 Octobev event had low counting rates and probably 
occurred partially outside the range of the OSO small raster. Therefore only 
the event occurring at U i 0259-0323 on 29 October was suitable for analysis. 

The event had just started at the OSO night-day transition, so the analysis 
covers the flare peak and decay phases. 

The Goddard data is output as one raster every 61.44 seconds. C, the 
difference in the total counts behind the two filters, is printed out. In addition, 
at each of the 368 readout points in the raster is recorded a "-2", a M 0 n or an 
even positive integer. If the manganese filter counts exceed the iron filter 
counts at "-2" is recorded. The even positive integers, n, represent a 
difference of n or n+1 counts with the iron filter count higher. Outside 
flaring regions, recordings higher than +2 are rare. Therefore readings of 
n = 4 or greater were regarded as real and counted as n+1/2 counts, if they 
occurred in the flaring region. Readings of 2 were counted as 2 counts in the 
flaring region. Readings of 0 and -2 were not counted. The UCSD counting 
rate varied slowly with time, and the Goddard rate was assumed to be constant 
during that part of the raster that covered the flaring region. 

To find the Goddard flux we first compute the sum of the individual counts 
c^ as outlined in the previous paragraph. We then need to subtract out back- 
ground at the n^ flare points in the raster. The total number of background 
counts at 368-n£ points is 

B = C- jJ c. , (4. 1) 

i=l 1 


so that the background at the n£ 



points is 



Then the total number of flare counts is 


(4.2) 
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Here T- is the temperature times 10 K, N Q is the electron density times 

-9-3 ' -19 

10 cm , I g y is (ground state ionization potential / 1.6 x 10” joules), £ - 1, 

and F(Z, z) = 1.8, Thus tj ~ 131 /N^ for T = 2 x 10^K. Similarly the recombina- 
tion time is given by 

t R = 10 4 T^ / ZN r/ (4.7) 

where Z = 26. At 2 x 10 K, t R = 544/Ng, Thus, since the continuum flux and 
temperature vary smoothly, 'he total flux in the Ross filter pass band should be 
a smooth function of time. We therefore made a least squares linear fit to 

= mt +b, with t measured in seconds from 0300 UT, and found m = -3. 865 and 
b = 8093. This curve is shown in Figure 5. The rms deviation of the data from 
this curve, 1384 photon-cm -sec” , is an estimate of the error in 
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The UCSD measurements were used to separate the line and continuum 
components of Oq. As noted above, UCSD channels 3 and 5-8 are strongly 
dominated by continuum. The Monte Carlo program waB used to find the best 
fit temperature and emission measure as described in Section 2. Using these 
parameters and the Aerospace model X-ray spectrum computer program, the 

O 

continuum flux in the 0, 174 - 0. 190 A range was computed and subtracted from 
to obtain the line flux, Input line counting rates of five and ten percent of 
the total UCSD integral rate were used. Comparison of the resulting line fluxes 
with those input to the UCSD program for consistency indicated that the lower 
line fluxes were more accurate. Since the continuum obtained in this way is 
essentially an interpolation of the UCSD data, we estimate the continuum sub- 
traction contributes at worst an error of 20% of the line flux. Table II shows 
the results of this computation, assuming a helium abundance of 0.20 to cal- 
culate the emission measure, and taking from the fit in Figure 5. The temp- 
erature was determined as an average of the values obtained for channels 5-8 
using plots such as Figure 3, assuming an isothermal plasma. The errors are 
the rms deviations from the quoted values. 

If a low density plasma is assumed, an estimate of the Fe XXV population, 
a 24 (T z ), and of T z may be obtained from the data. Values for Ap e and A^ e 
must be assumed. Combining equations (3.1), (3,2), (3.4), (3.7), and (3.9) 
gives for the total line flux. 


1 + G + ,2963F d (T) + .0808FJ 

-2 

photon-cm -Bee 

In equation (4.8), the factor converts from energy to photon flux under the 
approximation that all lines included have the 'tame wavelength. The ratio of 
the ion populations is introduced since the applicable one is a. A (T ) and ® (T) 
contains the factor Since me inner shell term is usually small, we begin 

by setting the line flux from Table II, equal to the remainder of equation (4, 8): 


(T z )] (4. 8) 

-1 J 


*T< T ’ T z> = 


* 3 < T > 


a 24< T z> 


E, 


R 


24 


(T) 



*s< T > 



*24< T z> 

a 24< T > 


[ 1 + G + .2963 F d 



(4.9) 


Table II 


Time 

*G 

, -2 -1. 
(cm *ec ) 

*c 

/ - 2 -1* 
(cm sec ) 

*L 

(cm'^.gec*^) 

T 

(10 6 K) 

EM 

(10*®cm‘ 

0301:45 

7687 

3279 

44u» 

21.2 ± 0. 9 

1. 31 

0302:53 

7424 

4543 

2881 

22.2 ± 1.0 

1.58 

0303:48 

7204 

4960 

2244 

22.2 ±0.7 

1 . 67 

0304:50 

6972 

4805 

2167 

21. 1 ± 1. 1 

1.79 

0305:57 

6713 

4394 

2319 

20.0 ± 0. 8 

2. 14 

0306:58 

6477 

4120 

2357 

19. 1 ± 0. 6 

2. 35 

0308:00 

6238 

3750 

2488 

17. 9 ± 0.3 

2.64 

0308:56 

6021 

3525 

2496 

17.4 ± 0.4 

2.87 

0309:57 

5786 

3158 

2628 

16. 5 ± 0, 3 

3. 18 

0311:04 

5526 

2553 

2973 

15. 1 ± 0.4 

3. 79 

0312:06 

5287 

2219 

3068 

14.4 ± 0.9 

4. 10 

0313;01 

5074 

2073 

3001 

14. 1 ± 0.9 

4. 26 

0314:09 

4812 

1775 

3037 

13. 3 ± 0. 8 

4.86 

0315:04 

4599 

1778 

2821 

13.5 ±0.9 

4. 50 

0316:11 

4340 

1378 

2962 

12.4 ± 1. 1 

5.47 

0317:13 

4100 

1233 

2867 

12.2 ± 1. 0 

5.40 

0318;08 

3888 

1106 

2782 

11.9 ± 0.9 

5.60 

0319:10 

3648 

923 

2725 

11.4 ± 0.5 

5.96 


; i 


With assumed values of a h. and A Fe and measured T and EM, equation (4.9) is 
solved for a 2 4 <T r )/ a 24 (T), Since a 24 (T) is known, a 24 (T z ) and, hence, T 
may be calculated. This value may then be used to find the inner shell excitation 
contribution to the flux and equation (4. 8) may be iterated to improve the estimates 

of a 24' T z^ and Ta- 
xable III summarizes the results of the calculation in the preceding paragraph 


for A Ffi = 2. 0 x 10 


-5 


and A^ e = 0. 2. 


For lower A^ e the emission measure in- 


creases as Figure 6, which plots the calculated continuum for A^ e = 0,2 and 
Ajj e = 0. 15 relative to that for A^ e = 0.07, illustrates. For higher A Fo' a 24< T *> 
is, of course, lower, since the line flux is proportional to Ap e X “24< T z>- 

Table Til also records the product T * EM. Assuming a constant flare density 
as a function of time, T* EM is a measure of the flare energy content since the 
thermal energy, E.J,, is 

E t = ~ <n e +Enj)VkT, (4.10) 

where V is the volume and the sum is over all ions present in the plasma. Since 


£n r » n H (l + A He ) , (4.11) 

and at coronal temperatures 

n e 53 n H (1 + * ( 4 * I2 ) 

we have 

E T » I n VkT ( 1 t ■■ \ , (4.13) 

2 \ 1+2A H ./ 

Finally we rewrite (4, 13): 

e t „i sm kT /i + l-..ta£\ . < 4 - 14 > 

2 n e \ 1 + 2A H J 

Table III indicates that the flare energy content continues to rise even as the flare 
cools. A similar result was reported by McKenzie, et al . (1973). Numerically, 
the peak energy is 2.6 x 10^/ n g ergs. 
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Table 111 


Time 

T 

(10 6 K) 

EM 

<10 4 Vffi' 3 ) 

EM* T 

<10 55 cm“ 3 K> 

*24< T > 

ft 24<V 

T z 

(10 6 K) 

0301:45 

21.2 

1.31 

2.78 

.693 

.501 

17.9 

0302:53 

22.2 

1.58 

3.51 

.726 

.21 D 

14.7 

0303:48 

22.2 

1.67 

3.71 

.726 

. 176 

14.0 

0304:50 

21.7 

1.79 

3.88 

.713 

. 150 

13.6 

0305*57 

20.0 

2. 14 

4.28 

.646 

. 168 

13.9 

0306:58 

19. 1 

2. 35 

4.49 

.594 

. 176 

14.0 

0308:00 

17.9 

2.64 

4. 73 

. 503 

.213 

14.4 

0308:56 

17.4 

2.87 

4.99 

.455 

.200 

14. 3 

0309:57 

16.5 

3. 10 

5.25 

.394 

.230 

14.7 

0311:04 

15. 1 

3. 79 

5.72 

.268 

.282 

15.2 

0312:06 

14.4 

4. 10 

5.90 

.209 

.314 

15. 5 

0313:01 

14. 1 

4.26 

6.01 

. 186 

.322 

15.6 

0314:09 

13.3 

4.8 6 

6.46 

. 125 

.336 

15. 7 

0315:04 

13.5 

4.50 

6.08 

. 140 

. 318 

15.5 

0316:11 

12.4 

5.47 

6. 78 

.076 

. 373 

16.0 

0317:13 

12.2 

5.40 

6.59 

.068 

.394 

16.5 

0318:08 

11.9 

5. 60 

6, 66 

.055 

.426 

17. 1 

0319:10 

11.4 

5. 96 

6. 79 

.036 

.405 

16. 7 


®g^ ,00 ‘ 


5. 


Result* 


The data were analyzed under two assumptions. First the electron density 
was assumed to be low enough that T z could differ from T. Second, an estimate 
of Ap e was made assuming = T. 

Table III in Section 4 presents some of the results obtained under the low 

density assumption. The time development of (T z ) does not occur as one would 

expect. For times at which T exceeds T^ the plasma should be ionizing; that 

is T should be increasing. Kafatos and Tucker (1972) show some calculated 
z 

ion population curves for hypothetical X-ray flares. Because of the large probable 
error in 4^., the data in Table III may not be inconsistent with the expected 
theoretical flare decay. 


A second approach was taken to comparing theoretical low density cal- 
culations to the experimental data. For any ion, the ionization balance equation 
is (Kafatos and Tucker, 1972) 


da z 

■ar = n e fa .+i Vh 


a a, + Q .a ,-Qa) 
z z z-1 z-1 z z 


(5. 1) 


In equation 5. 1, a z is the fractional population of the ion with charge z, a z is the 

total recombination coefficient for the process z -* z-1, and Q z is the ionization 

coefficient for z -* z+1. Iron near or in a solar flare may be in any state from 

~ Fe XVI to Fe XXVI so that to obtain a is a function of time a large number of 

z 

coupled equations must be solved. We have made a number of simplifying 
assumptions in treating equations (5. 1). We assumed that only three stages 
of ionization are important: Fe XXIII, Fe XXIV, and Fe XXV. All iron in 
states below Fe XXIII was added to that in Fe XXIII, This is not too bad an 
approximation since the ionization time for Fe XXII is only ~ 50/N^ seconds. 

The approximation slightly increases the Fe XXV ionization rate. The second 
approximation was to ignore recombination from Fe XXVI to Fe XXV. At the 
temperatures, T, of the flare in question this is an excellent approximation 
because the population of Fe XXVI is very small. As a result of the approxi- 
mations we were left with three coupled linear differential equations in three 
unknowns. These were solvable by matrix methods without too much labor. 

The a' s and Q’ s were obtained from Jordan (1970). 
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Equations (5. 1) were solved at each readout interval from 0304:50 on. 

At w*ach time, the increase in emission measure from the previous readout was 
considered to be due to an increase in flare volume at constant density. All 
iron included in this increase was considered to be in a low ionization state 
and was deleted when computing the input to the next calculation. Thus if EM 
increased by 5%, the values of a^, a and obtained from the nth solution 
of equations {5. 1) were multiplied by (1.05)“* and used aa initial conditions for 
the (n+l)th solution cf the equations. The initial conditions at 0304:50 were set 
at a^j 3 .641, a^j ” .129, and a^^ = .161. The electron density n fi was set at 
10 9 , and Ap e = 2 x 10“^ and A^ e = 0.2 were used. 

The results of the foregoing calculation are shown in Table IV, T^ is 
given as the temperature at which a£4(t) would obtain if the plasma were in a 
state of ionization equilibrium, 4^ (T, T z ) is 4 rp from equation (4.8), and 
4 is the continuum in the 0. 174-0, 190 band discussed in the previous sections. 
Thus 4^ + 4^ should be compared with the raw data 4^. This comparison is 
shown in figure 7. The agreement between the data and the calculation appears 
to be good except that the Goddard flux persists after the calculated curve has 
fallen off. The UCSD data were not available for this period, but it is likely 
that they would be unusable, since as T decreases the proportional counter data 
become unreliable. 


Despite the apparent agreement between data and calculation in figure 7 

the low density solution is probably not valid. From the rasters an upper limit 

20 2 

on the flare area of about 10 cm may be set. The ATM data for a large 

variety of flares suggest that the actual area is much smaller (Vorpahl et al. , 

1975; Kahler et al. , 1975 ). With EM » 3 x 1 0^** cm \ we have n h ^ 3 x 10^, 

9 e 

Since flare heights, h, are typically ~ 10 cm (Catalano and Van Allen, 1973), 

2 19 q 

n «* 3 x 10 and n ^ 5 x 10, The solution to equations (5. 1) indicates that 

e e 9 

at such high densities, T is much nearer T than it is for n = 10 . Therefore 

z e 

the high density case (T = T) must be examined. 

The high density case is easy to treat. One simply solves equation (4.8), 
assuming T = T , To determine the iron abundance, we calculate $ w » defined 
as 4^, in equation (4. 8) for Ap e = 10"^. We then make a least squares fit of 
the form m4 = 4, where 4 T is the measured line flux, 4 (measured) - 4 : 

W JLj V*/ C 

that is we minimize 
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Table IV 



4 


J 



Time 

*24< l > 

T z at 
a 24 (t) 

* L (T f T z ) 

* 

c 

* + $ 
Ld 

aE 

1 

0304:50 

. 161 

13.8 

2070 

4805 

6875 

i 

t 

0305:57 

.220 

14.5 

2710 

4394 

7104 

1 

0306:58 

.268 

15.1 

3209 

4120 

7329 

f 

1 

f. 

0308:00 

. 305 

15.4 

3410 

3750 

7160 

1 

i: 

j 

It 

0308:56 

. 306 

15.4 

3470 

3525 

6995 

0309:57 

.278 

15.2 

2948 

3158 

6106 

r ; 

t 

0311:04 

.289 

15. 3 

2885 

2553 

5438 

U 

i» • 

I 

j* r * 

0312:06 

. */5 

15.2 

2569 

2219 

4788 

k 

i 

0313:01 

.265 

15. 1 

2372 

2073 

4445 

i. 

0314:09 

. 234 

14.7 

2046 

1775 

3821 

j; 

0315:04 

.247 

14.9 

2103 

1778 

3881 

* 

0316:11 

. 203 

14.3 

1597 

1378 

2975 

i' 

0317:11 

. 195 

14. 2 

1397 

1233 

2630 


0318:08 

. 155 

13.7 

1027 

1106 

2133 




19 - 


(5.2) 


=X>*»i - * L1 > : 


Differentiating, we obtain 


m 


? 


*Li *wi 


5*^ 

Assuming that the error in 4^ dominates we estimate the error, 6m, as 


(5.3) 


6m = |T^|(|HL\ 2 (64) 2 t . 1 


or 


6m 


im 


]_ 

2 


(5.4) 


2 


Lj 


(5.5) 


wj 


From equations (5, 3) and (5. 5) we obtain m = 1,11 £ . 14, for A He = 0. 2. For 
lower values of A^ e , m will be lower in inverse proportion to the assumed emission 
measure. The assumed emission measure is in turn inversely proportional to 
the curves in figure 6, Thus for Ajj e = 0,07, m« 0.8. For Apj e = 0,2 we 
estimate Aj, e = 1. 1 x 10 " 5 . The data are displayed fn Table V. 

Figure 8 shows a plot of 4 + 1.11 4 compared to the data, 4_. The 

C W Lt 

point at t = 173 was excluded from the calculation of m because of its great 
deviation from the points around it. As in the low density case the curve predicts 
less iron line flux than is observed late in the flare. 
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Table V 


Time 

*wl 

*Gi 

*ci 

*Li 

0301:45 

2878 

8394 

3279 

5100 

(0302:53) 

(3950) 

(4214) 

(4543) 

(-329) 

0303:48 

4175 

8994 

4960 

4034 

0304:50 

4230 

8055 

4805 

3250 

0305:57 

3772 

7555 

4394 

3161 

0306:58 

3427 

8824 

4120 

4704 

0308:00 

2747 

5417 

3750 

1667 

0308:56 

2532 

4386 

3525 

861 

0309:57 

2061 

7063 

3158 

3905 

0311:04 

1340 

5907 

2553 

3354 

0312:06 

985 

6472 

2219 

4253 

0313*01 

842 

2496 

2073 

423 

0314:09 

541 

5008 

1775 

3233 

0315:04 

606 

4249 

1778 

2471 

0316:11 

308 

2807 

1378 

1429 

0317:13 

250 

4122 

1233 

2889 

0318*08 

187 

3585 

1106 

2479 

0319:10 

126 

1711 

923 

788 
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6. Discussion 


-5 

In the previous section we arrived at an iron abundance of I, I x 10 for a 

helium abundance of 2 x 10~* or an iron abundance of ~8 x 10"^ in the limit of 

very low helium abundance. These values lie considerably below those found 

in other investigations. Walker et al. (1974a) recently determined an iron abun- 
.5 ■» — 

dance 2.6 x 10 for a solar active region, based upon a detailed analysis of 

Fe XVII emission. Walker et al. (1974) had earlier placed the iron abundance at 

4 x 10- 5 , based upon a less complete model of the Fe XVII ion and lower quality 

observations than were used in the later paper. On the low side. Dupree (1972) 

-5 

se«. the iron abundance at 2. 0 x 10 based on EUV observations. 

Can our observations be reconciled with previous iron abundance estimates? 

To answer this question we must consider possible systematic errors and un- 
certainties in the calculation of the expected iron line flux. As we note in Section 
III, we have probably overestimated a, the parameter relating the unresolved 
dielectronic recombination satellite line flux to the resonance line flux. If so. 
the expected line flux for a given Ap Q iB overestimated by a factor of about 1. 1. 

We noted in Section IV that the continuum subtraction might result in an error 
as large as 20% in the line flux. Thirdly a systematic error in the temperature 
estimate in the early phase of the flare could result in an error in the expected 
flu. 1 ';. For example, an overestimate of T by 2 x 10^K at a true temperature of 
2 x 10 K results in an underestimate of EM by 17% and an overestimate of the 
expected line flux per unit EM by 51%, for a net overestimate of the expected 
flux by 29%. Since the high counting rates are heavily weighted, a temperature 
error here would be most effective. Errors larger than that discussed above 
aie improbable. If all of the above factors are considered, the result for A^ e = 

0. 2 becomes (1. 1 t . 14) x (1. 1) x (1.2) x (1.29) = 1.9 t .24. 

Two other factors are not as easy to estimate. Sufficient high quality iron 
line spectra are not available to verify theoretical line strength calculations. The 
factor, G, is based upon measurements of low Z ions, Parkinson (1972), however, 
from an excellent magnesium spectrum, found G = 1.04, in good agreement the 
oxygen measurements of Rugge and Walker (1971) and with our assumed value of 

1. 1. The satellite line intensities have not been verified, either. Finally, we 


V 
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should consider the consequences of a breakdown of the assumption that the 
flare plasma is isothermal. Because the continuum flux increases more than 
in proportion to the temperature, if the plasma contains components at temp- 
eratures greater than T (as it must if the plasma is not isothermal), the em- 
ission measure based upon an isothermal plasma will be an overestimate. 

7 

However, the line flux increases more rapidly with T (around 2 x 10 K) than 
does the continuum. Then the line to continuum ratio will be higher than that 
for an isothermal plasma. Thus a breakdown in the isothermality assumption 
would result in an overestimate of the iron abundance. 

Taking into account the discussion of the preceding paragraphs, it appears 

that our may be reconciled with other low estimates but not with values 
* ® -5 

much in excess of 2 x 10 . Further, the helium abundance must be high for 

this reconciliation, This experiment thus lends support to a helium abundance 
relative to hydrogen around 0.2. The measurement is not good enough to make 
a definitive determination of Apj e relative to Ap e< 

The determination of the solar helium abundance is an important unsolved 
problem. The continuum emission processes are well enough understood that 
the helium abundance is the major unknown in calculating continuum emission 
in some wavelength and temperature regimes (see Culhane, 1969). We have 
chosen to attack the problem by measuring the line to continuum ratio for 
helium -like iron during a solar flare, in the hope that a definitive iron abundance 
measurement, coupled with our observations, would allow the helium abundance 
to be determined. This approach has the disadvantage that iron line emission is 
not well understood and that solar flare conditions impose analysis difficulties. 
However, flare observations have the advantage that substantial continuum emis- 
sion occurs in a short wavelength region that is line free. A large number of 
high quality solar flare iron line spectra need to be analyzed to verify the iron 
line emission theory and the iron abundance before the technique used here can 
yield a helium abundance determination. Crystal spectrometer or Ross filter 
measurements, together with high quality continuum spectra, perhaps measured 
with a cooled solid state detector, can then provide a good determination of the 
solar helium abundance. 
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FIGURES 


Figure 1 
Fir't e 2 

Figure 3 

Figure 4 
Figure 5 
Figure 6 

* xgure 7 
Figure 8 


The man absorption coefficients of iron and manganesr. 

The calculated response of the UC5D proportional counter 
to thermal continuum spectra. The effective channel edges 
reflect a 6% gain decrease »n this plot. 

The UCSD detector response at a total counting rate of 
4-1 3-1 

2 X 10 sv*c , J X 10 sec - of which are iron line 
cjunts. A 5% gain decrease and an additional 5% decrease 
due to high counting rates are assumed. 

The temperature dependent multiplying factors for inner- 
shell excitation and dielectronic recombination sa ijllite 
lines. Data are from calculations by Bhalla et al. (1975). 

The Goddard 0.174 - 0. 190 nm x-ray flux as a function 
of time and the linear least squares fit to the data. One 
±ltr error bar is shown. 

The continuum at 0.2 nm as a function of electron tempera- 
ture with the helium abundance. A. ? , as a parameter. The 

flux is normalized such tnat the flux for = 0. 07 is unity 
at each temperature, 

Comparison of 4. + 4> ( calculated assuming density = 

a _3 " c 

10 cm with the measured 0. 174 - 0. 190 nm flux. 

Comparison of * +1.11* and for the case T = T . 

The helium abundance is assumed to be 0.2, relative to 
hydrogen. 
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